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ABSTRACT 

The present paper is concerned with the application of an existing 
analytical crack-closure model to study crack growth under various load his- 
tories. The model was based on a concept like the Dagdale modal, but modified 
to leave plastically-deformed material in the wake of the advancing crack tip. 

The model was used to correlate crack-growth rates under constant-amplitude 
loading, and to pr.^dict crack growth under variable-umplitade and aircraft- 
spectrum loading on 2219-T851 aluminum alloy sheet material. The predicted 
crack-growth lives agreed well with experimental data. For 80 crack-growth 
tests subjected to various load histories, the ratio of predicted-to-experimental 
lives (Np/N,p) ranged from 0.5 to 1.8. The mean value of Np/N^ was 0.97 and 
the standard deviation was 0.27. 
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material crack-growth constants (k ■ 1,5) 
half-length of crack, m 
half-length of final crack, m 
half-length of initial crack, m 
half-length of starter notch, m 

half-length of crack plus tensile plastic zone, m 

boundary-correction factor on stress intensity 

1/2 

maximum stress-intensity factor, MPa - m 

1/2 

elastic-plastic fracture toughness, MPa - m 
fracture toughness parameter 
number of cycles 

number of crack-growth delay cycles 

number of cycles predicted from analysis 

number of cycles from test specimen 

stress ratio 

applied stress, MPa 

maximum applied stress, MPa 

minimum applied stress, MPa 

crack-opening stress, MPa 

specimen thickness, m 

specimen width, m 

constraint factor, rt = 1 for plane stress 
and ot - 3 for plane strain 
stress-intensitv factor range, MPa - m^/^ 
effective stress-intensity factor range, MPa - m 


1/2 


1/2 

ef£ectlv« threshold atress'^lntenslty factor range, MPe - m. ' 
AK^j^ threahold atreaa~lntenalty factor range, MPa - 
ASgff effective stress range, MPa 

p length of tensile plastic lone, m 

plastic-aone siae calculated from overload, m 
flow stress (average between and o^), MPa 

Oy^ yield stress (0.2 percent offset), MPa 

0^ ultimate tensile strength, MPa 

u length of compressive plastic zone, m 


INTRODUCTION 

Fatigue cracks remain closed during part of the load cycle under fatigue 
loading. The crack-closure concept has been used to correlate crack-growth 
rates under constant-amplitude loading tl»2] and is a significant factor in 
causing load-interaction effects on crack-growth rates (retardation and 
acceleration) under variable-amplitude loading. Fatigue-crack closure is 
caused by residual plastic deformations remaining in the wake of an advancing 
crack. 

The crack-closure phenomenon has been analyzed using two-dimensional, 
elastic-plastic, finite-element methods [3-6). The finite-element analyses 
were shown to be quite accurate, but were very complicated and required large 
computing facilities. There nave also been several attempts to develop simple 
analytical models of crack closure [3,7-12). All of these r >dels were based on 
a concept similar to the Dugdale model [13) or strip-yield model, but modified 
to leave plastically-deformed material in the wake of the crack. Newman [3j, 
Budiansky and Hutchinson [8), and Fuhring and Seeger [10,11) studied only the 



crack-closure behavior. But, Dill and Saff [7], Hardrath* Navaan, Elber and 
Poe [9], and Newman [12] used the crack-opening stresses from the models to 
predict crack growth under spectruip loading. 

The purpose of the present paper Is to apply an existing analytical crack- 
closure model [12], which simulates ^ lane-stress and plane-strain conditions, 
to crack growth under various load histories. The model was based on the 
Dugdale model [13], but modified to leave plastically-deformed material along 
the crack surfaces as the crack advances. Plane-stress and plane-strain condi- 
tions were simulated by using a "constraint" factor on tensile yielding. 

The crack-closure model, developed in reference 12, was for a central 
crack in a finite-width plate that was subjected to a uniformly applied stress. 
To calculate Elber 's effective stress-intensity factor range [2], crack-opening 
stresses were calculated from the model under constant-amplitude loading at 
various applied stress levels and stress ratios. Experimental crack-growth rate 
data from d219-T851 aluminum alloy sheet material under constant-amplitude 
loading [14] were correlated with the effective stress-intensity factor range 
for a wide range of stress levels and stress ratios. An equation relating 
crack-growth rate to effective stress- intensity factor range, threshold stress- 
intensity factor range, and fracture toughness, developed in reference 12, was 
applied herein over the total range of crack-growth rates. The closure model 
was then used to predict crack growth in 2219-T851 aluminum alloy sheet 
material under variable-amplitude and aircraft-spectrum loading [14]. 

ANALYTICAL CRACK-CLOSURE MODEL 

The following section is a brief description of the analytical crack- 
closure model developed in reference 12. 
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To calculate crack-closure and crack-opening stresses during crack propaga- 
tion, the elastic-plastic solution for stresses and displacenents in a cracked 
bod)? must be knovm. Because there are no closed-form solutions to elastic- 
plastic cracked bodies, simple approximations must be used. The Dugdale model 
[13] is one such approximation. The crack-surface displacements, which are used 
to calculate contact (or closure) stresses under cyclic loading, are influenced 
by plastic yielding at the crack tip and residual deformations left in the wake 
of the advancing crack. The applied stress level at which the crack surfaces 
become fully open (no surface contact) is directly related to contact stresses. 
This stress is called the "crack-opening stress." 

The model was developed for a central crack in a finite-width specimen sub- 
jected to uniform applied stress, as shown in figure 1. The model was based on 
the Dugdale model but modified to leave plastically-deformed material in the 
wake of the crack. The primary advantage in using this model is that the 
plastic-zone size and crack-surface displacements are obtained * superposition 
of two elastic problems; a crack in a finite-width plate subjected to (1) a 
remote uniform stress, S, or (2) a uniform stress, O, applied over a segment 
of the crack surface. The stress-intensity factor and crack-surface displace- 
ment equations for these loading conditions are given in reference 12. 

Figure 2 shows a schematic of the model at maximum and minimum applied 
stresses. The model was composed of three regions: (1) a linear elastic region 

containing a fictitious crack of half-length c + p, (2) a plastic region of 
length p, and (J) a resldvial plastic deformation region along the crack sur- 
faces. The physical crack is of half-length c. The compressive plastic zone 
Is <jj. Region 1 was treated as an elastic continuum, and the crack-surface 
displacements under vari^.is loading conditions are given in reference 12. 
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Regions 2 and 3 were composed of rigid-perfectly plastic (constant stress) bar 
elements with a flow stress » which is the average between the yield 

stress, Oyg, and the ultimate tensile strength, 0^. The shaded regions 
in tlgure 2(a) and 2(b) Indicate material which is in a plastic state. At any 
tpplied striiss level, the bar elements are either intact (in the plastic zone) 
or broken (residual plastic deformation). The broken elements carry compressive 
loads only, and then only if they are in contact. The elements in contact yield 
in compression when the contact stress reaches “Og* Those elements that are 
not in contact do not effect the calculation of crack-surface displacements. 

To account for the effects of state-of-stress on plastic-zone size a constraint 
factor a was used to elevate the tensile flow stress for the intact elements 
in the plastic zone. The effective flow stress otOg under simulated plane- 
stress conditions was and under simulated plane- train conditions 

was 30 q. The constraint factor is a lower bound for plane stress and an 
approximate upper hound for plane strain. These constraint factors were veri- 
fied using ela.itlc-plastic finite-element analyses of cracked bodies under plane 
stress [6] and plane-strain conditions. The procedure used to establish the 
constraint factor (a) used herein is discussed later. 

The analytical crack-closure model, discussed in detail in reference 12, 
was used to calculate crack-opening stresses, S^, as a function of crack length 
and load history. In turn, the crack-opening stress was used to calculate the 
effective stress -intensity factor range, as proposed by Elber, and, consequently, 
the crack-growth rates. 

FATIGUE-CRACK GROWTH RATE EQUATION 

The crack-growth equation proposed by Elber (2] states that the crack- 
growth rate is a power function of the effective stress-intensity factor range 
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only. Later, Hardrath, Newman, Elber, and Poe [9] shoved that the power law 
was Inadequate at high growth rates approaching fracture. The results presented 
in reference 12 showed that it was also Inadequate at low growth rates 
approaching threshold. To account for these effects, the power law was modified 
In reference 12 to 


1 - 


^ “ Cl AKgff 
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AK 
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where 


AKq = 1 + 


max 


( 2 ) 


and 


~ ^max ^ 


AK 


.ff ■ - So'; S 


(3) 


(4) 


The crack-opening stresses, S^, were calculated from the analytical closure 
model. Equation (1) gives the "sigmoidal" shape commonly observed when fatigue 
crack-growth rate data is plotted against stress-intensity factor range. In 
the intermediate range of crack-growth rates, equation (1) is basically Elber* s 

i C2 

I j proposed power law, Ci AK^^^. The constants Ci to C 5 were determined to 

I best fit experimental data under constant-amplitude loading. 

The coefficients C 3 and were determined from threshold data on the 

2219-T851 aluminum alloy sheet material from reference [15]. The effective 
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threshold stress'lntenslty factor range, was determined from the threshold 

stress-intensity factor range, ^th* 

“ “th ■ 1 . ■ r* “‘th (5) 

The coefficient C 5 is the elastic stress-intensity factor at failure or 

cyclic fracture toughness. The coefficient C 5 was chosen to be 

1/9 2/2 

77 Mia - m (70 ksi - in ' ) on the basis of the crack-growth tests in 
reference [14], (See ref. 12.) 

The coefficients and C 2 were found from constant-amplitude rate 

data [14], after C^, C^, and C 5 were determined, by using a least-squares 

regression analysis. The constant-amplitude correlations were made using Sq 
values computed from the model with various constraint factors. It was found 
that an a of about 1.9 would give a good correlation under constant-amplitude 
loading. The procedure used to obtain a will be discussed later. A summary 
of the coef ficiu..is used to correlate the constant-amplitude data with a ■ 1.9 
are as follows: 

= 2.486 X 10"^^ C- 314 x 10"®) 

C 2 = 3.115 

Cj = 2.97 MPa - m^^^ (2.7 ksi - in^''^) 

= 0.8 

C 5 =• 77 MTa - m^^^ (70 ksi - In^^^) 
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When SI units are used, ^eff Snax ^ dc/dN 

is given in a/ cycle. When U.S. Cuatomary units are used, and 

are given in ksi - in^'^^ and dc/dN is given in in. /cycle. 

Figure 3 shows a plot of AK against dc/dN for several R ratios for 
2219-T851 aluminum alloy sheet material to Illustrate the sigmoidal shape of 
equation (1). The experimental data were obtained from reference [15] and the 
curves were calculated from equation (1). The R * -1 data were obtained 
from a small center-crack tension specimen (W = 76.2 mm) and the other data were 
obtained from small compact specimens (W ■ 50.8 mm). The crack-growth coef- 
ficients (Cj^, Co, C 3 , and C^) used to calculate the curves were identical to 
those shown in equations ( 6 ). However, the coefficient C 5 for the small com- 

1/2 

pact specimens was 38,5 MPa - m ' and for the small center-crack specimen 

1/T 

was 55 MPa - m The coefficients (C 5 ) were calculated from the Two-Parameter 

Fracture Criterion [16] using Kp ■= 550 MPa - and ra * 1. These values 

of 5(p ana m were obtained from the final crack lengths and maximum stress 
levels used in the constant-amplitude tests from reference [lA], (See ref. 12.) 
The crack-growth rate curves in figure 3 are in good agreement with the experi- 
mental data. 

APPLICATION OF THE CRACK-CLOSURE MODEL AND RATE EQUATION 

The analytical crack-closure model [12] and crack-growth program (FAST - 
Fatigue Crack Growth Malysis of Structures) was applied to constant-amplitude, 
variable-amplitude and aircraft-spectrum loading on 2219-T851 alum^ m alloy 
sheet material. 

Under constant-amplitude loading, the model was exercised under simulated 
plane stress, plane strain, and ccx:ditions between these limits. The particular 
constraint factor (a) used herein, to approximate the state-of-stress, was 
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obtained from the constant-amplitude crack-growth rate data. The crack-growth 
rate equation (eq. (1)) was also determined from the constant -amplitude data. 

The sane constraint factor was also used to predict crack growth under 
variable-amplitude and aircraft-spectrum loading. The crack-opening strtj 
were calculated from the model as a function of crack length and load hiscory, 
and the crack-growth rates were predicted from equation (1). The predicted 
crack-growth lives are compared with experimental data in the following sections. 

CONSTANT-AMPLITUDE LOADING 

Crack-opening stresses .- Reference 12 showed that the calculated crack- 
opening stresses under constant-amplitude loading were independent of the con- 
straint factor for stress rati s (R) greater than about 0.7 and were equal to 
the minimum applied stress. Thus, AKgff is equal to AK for R > 0.7. 

Using crack-growth rate data from references 14 and 15 for R > 0.7, the crack- 
growth constants C^ and C 2 were determined to best fit the high R value 
data only. Basically, this crack-growth rate equation (eq. (1)) depicts the 
relation between snd crack-growth rate. If ^gff dc/dN are 

unique, then the crack-growth rates for tests at R ratio less than 0.7 should 
indicate the experimental and, consequently, the experimenta] value 

of Sq. This value of Sq is referred to as the "semi-empirical" r rack-opening 
stress. For each test, le semi-empirical Sq value was .. 'turned to be constant 
and was determined from a least-squares regression analysis. These crack- 
opening stresses, normalized by the maximum applied stress, are shown in fig- 
ure 4 as a function of the R ratio (symbols). The open symbols and the solid 
circular symbol are results from center-crack tension (CCT) specimens [14,15]. 

The solid triangular symbols are results from compact specimens [13]. These 
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results indicate values that correlate crack-growth rate data at various 
R ratios with the results at R > 0.7. The dashed line indicates where Sg is 
equal to or where is equal to AK. 

For R ratios less than 0.7» the calculated crack-opening stresses from 
the closure model are a function of the constraint factor. The semi-empirical 
results at R » 0 were used to estimate the constraint factor. A value of 
1.9 was found to give good agreement between the calculated and seml-emplrical 
values. The solid curves show calculations from the closure model with a ■ 1.9 
for various applied stress levels. The overall agreement between the calcu- 
lated and semi-empirical values was considered reasonable. 

Crack-growth calculations .- The crack-closure model with a •= 1.9 was 
used to calculate the crack-growth lives from the initial crack length Cj_ to 
the final crack length c^ for the constant-amplitude tests [14] used in 
figure 4. The ratio of predicted-to-experimental lives (Np/N ranged from 0.6 
to 1.8. The mean value of Np/N<j. ras 1.01 and the standard deviation was 0.32. 
These calculations were considered reasonable, in view of the scatter that 
occurs in fatigue-crack growth rate tests. 

VARIABLE-AMPLITUDE LOADING 

Crack-opening stresses .- The closure model (a = 1.9) was used to study the 
crack-opening stresses under various load histories. The calculated crack- 
opening stresses under single-spike and two-level loading are shown in figure 5 
as a function of crack length. Under spikr loading, Sq takes a sudden drop 
when the crack-tip region blunts due to the spike loading. As the crack grows 
into the overload plastic zone (Pql)» the values rapidly Increase until 

they reach a maximum value at about one-half of Pcii^ This is the point of 
minimum crack-growth rate. The Sq values then drop and approach the 
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stabilized crack-opening stress (dashed line) for the low-level constant- 
amplitude loading. The retardation effects (Sq greater than dashed line) are 
nearly eliminated when the crack has grown about one overload plastic-zone 
size. In contrast, under two-level loading the high load was applied from the 
initial crack length (c^ 3 mm) for about 2500 cycles. Again, the values 

rapidly increase as the crack grows into the overload plastic zone, but they 
reach higher values than those occurring onder the single-spike loading. Thus, 
retardation effects are much stronger after multi-overloads than after single- 
spike loading. 

The calculated crack-opening stresses under compression-tension and 
tension-compression spike loading are shown in figure 6 as a function of crack 
length. The horizontal line is the crack-opening stress for the high R ratio 
(0.5) constant-amplitude loading. Under compression-tension spike loading, 
the "compressive underload" (single downward load excursion) had no influence 
on the subsequent values. Bul the tensile overload caused the Sq values 

to drop immediately, then rapidly rise above those from the steady 
state constant-amplitude loading (dashed line) for about one overload plastic- 
zone size. If the compressive underload is applied immediately after the 
tensil.' verload, the Sq values are considerably lower than those from the 
compression-tension spike loading. Thus, the compressive underload after the 
spike eliminates some of the retardation effects due to the tensile overload. 

A larger compressive underload after the spike causes a larger reduction in 
the values, but did not completely eliminate the retardation effect due 

to the tensile overload. 

Figure 7 shows the calculated crack-opening stresses during repeated 
compression-tension spike loading as a function of crack 1 .gth. The first 
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load level (1 >■ 1) was applied for 2500 cycles. The compression-tension spike 
(i * 2) was then applied. These load sequences were repeated until the speci- 
men failed. The dashed line shows the stabilised crack-opening stress 
for level 1 only. The dash-dot line shows the stabilized crack-opening stress 
if level 2 only was applied. These results show that the interaction bet%>een 
levels 1 and 2 cause values to increase slightly during the compressive 

underload, drop abruptly during the tensile overload and rapidly increase 
during the application of the 2500 cycles. Again, the values reach a 

maximum as the crack grows into the plastic zone caused by the tensile over- 
load. The Sq values would approach if the tensile overload was not 

repeated. 

Crack-growth predictions .- The crack-growth rate, at each load cycle, was 
computed from equation (1), using the current values of and S^. 

Equation (1) predicts retardation (or acceleration) if is larger (or 

smaller) than the crack-opening stress that would have been produced under 
constant-amplitude loading at and S ^.,- „ . To demonstrate how crack- 

grouth rates were calculated under variable-amplitude loading, an example is 
given. Figure 8 shows a typical variable-amplitude load histroy. The growth 
rate was computed from equation (1) using 

^eff " ^ 


wnere effective stress range on the k*"^ cycle. The growth 

increment per cycle is 


Ac 


k " ( 


dc 

dN 


( 8 ) 
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On the first and second tensile load excursion* S, 


®ink 


^^effk " ^naxk ~ 


where k >- 1 or 2, respectively. This equation was proposed by Elber [2]. 
However, on the third tensile load excursion, Is greater than 5^^, 

therefore the effective stress range was assumed to be 

1 


^Seff, 


(^maxk “ ^o) “ (^mink ~ ^o) 


C 2 " 


( 10 ) 


where k “ 3 and C 2 is the power on the growth law. Thus, the growth 
increment, Ac 2 + Ac^, is slightly larger than Ac]^, if The 

use of equation ( 10 ) was necessary because no crack-growth law, when expressed 
in terms of a power function (C 2 i 1 ), would sum to the correct growth incre- 
ment under variable-amplitude loading. For instance, if the load excursion 
^max 2 ^min^ extremely small, then the sum of growth Increments Ac 2 

and Ac 3 should be equal to the growth increment Acj^. If was less 

than Sq, then the growth increment Ar .3 •should be equal to the growth incre- 
ment Ac^. Equation (10) accounts for these limiting behaviors. Equation (10) 
is applied only when is greater than and only when the current 

maximum applied stress is higher than the highest maximum stress occurring 
since a stress excursion crossed S^. On the fourth excursion, ^Sg££ was, 
again, corouted from equation (9). The effective stress range on the 5th 
and 7th excursion were, again, computed from equation (10). But on the 6 th, 
8 th, and 9th excursion. 


° ^inaxj^ " ^mink 


( 11 ) 
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Equation (11) 


where k ■ 6, 8, or 9, respectively. Note that ^max^ - ®max 5 * 
was also proposed by Elber [2], 

Figure 9 shows the effect of the number of overload cycles on predicted 
crack-growth delay. Crack-growth delay Is the additional number of cycles 
required to grow the crack to fallurst following an overload, over the nuod>er 
required to grow the crack to failure under constant-amplitude loading only 
(level 3). The predicted results (symbols) show that delay or retardation is 
longer for larger nu]id>er of overload cycles. These results are in quantitative 
agreement with experimental observations [17] . Other retardation models 
account for the effect of the number of overload cycles on crack-growth retarda- 
tion empirically [18], or do not account for it at all [19,20], 

A comparison between experimental and predicted crack-length-against- 
cycles curves during two-level loading is shown in figure 10. The load sequence 
is shown in the insert. The high load (level 1) was a factor-of-2 larger than 
level 2. The predicted results (solid curve) were calculated from equation (1) 
using the values computed from the closure model. The maximum computed 

value of S^, during the application of level 2, was about 115 MPa and the 
minimum crack-growth rate was about 1.4 x 10”^ m/cycle. The predicted life 
was about one-half of the experimental life (symbols). The dash-dot curve 
shows the predicted results using no load interaction. 

Figure 11 shows a comparison of experimental and predicted crack-length- 
against-cycles curves for repeated tension-compression and compression-tension 
spike loading. The load sequences applied are shown in the inserts. The 
experimental results (symbols) and the predicted results (solid curves) show 
that the compressive underload applied after the tensile overload causes the 
crack to grow faster than when the compressive underload occurs before the 
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tensile overload. Although the predicted results show a stronger effect of 
the cooqpresslve underload than the axperliaental data, the agreeaents between 
the predicted and experiisental data are considered good. 

A conparison of experiaental (symbols) and predicted (solid curve) crack- 
length-agalnst-cycles curves for a repeated block loading sequence is shorn in 
figure 12. The load sequence is shown in the Insert. In contrast to the 
previous case, fifty cycles, of the tensile overload were applied before the 
compressive underload. The predicted results are in good agreement with the 
experimental data. 


SPECTRUM LOADING 

Crack-opening stresses .- The variation of crack-opening stress with crack 
length for a typical spectrum loading test is shown in figure 13. The half- 
length of the elox notch (cjj) was 3.2 mm. The specimen was cycled under 
constant-amplitude loading (S^^x “ MPa) at R » 0 until the crack grew to 
a crack half-length (c^) of 3.8 mm. Next, a typical fighter aircraft spectrum 
was applied to the specimen. The maximum stress was about 183 MPa and the 
minimum stress was about -30 MPa. The particular spectrum loads applied are 
given in refereni. [14] under test M91. The calculated crack-opening stresses 
plotted In figure 13 show only a small fraction of the number of values computed 
from the model. The crack-opening stresses follow a very irregular pattern 
while the cyclic loads are applied; even so, they tend to oscillate about a 
mean value. 

The Use of an "equivalent" crack-opening stress concept would greatly 
reduce the computer times required to complete a simulated test. The use cf 
an equivalent stress is justified because, at low to medium stress levels, the 
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crack-opening stresses stabilize under constant-anplitude loading* They also 
tend to oscillate about a nean value under spectrun loading. The equation for 
:he equivalent crack-opening stress, Sq, was 


E Ack" 


( 12 ) 


where the sunanation was perforated over the crack extension increaient 
Cf + 5 to + 10 Pjpax* maximum plastic-zone size, Pn»ax» 

calculated using the maximum stress in the spectrum. (For extremely high stress 
levels and low R ratios, where values do not stabilize, the simulated 

test specimen may fail before the equivalent crack-opening stress routine is 
activateu.) The dashed line in figure 13 shows the calculated equivalent 
crack-opening stress. The predicted crack-growth life using Sq was 3.5 per- 
cent less than the predicted life using S^, but the computer time was only 
about one-half as large (2.6 minutes to 5.6 minutes). 

C* ck-growth predictions .- In reference 14 crack-growth tests were con- 
ducted on center-crack tension specimens subjected to five basic aircraft-type 
load spectra. Three of the spectra were each applied at three different scale 
factors (same shape spectrum with different scaling of the stresses), and the 
other tv ^ spectra were each applied at two different scale factors. There 
wc-"' -hirteen different spectrum loading tests. 

Figure 14 compares predicted and experimental crack-length-against-cycles 
curves for a typical fighter spectrum. The specimens were subjected to the 
same spectrum, but with three different scale factors (0.2, 0.3, and 0.4). The 
predicted results using a = 1.9 (solid curves) are in good agreement with the 
experimental data (symbols). However, for all spectrum tests conducted at a 
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low-stress level (scale factor ■ 0.2), the predicted results gave longer lives 
than the experimental data. At the low-stress level, the plastic-sone slses 
are small compared to thickness and plane-strain conditions may prevail. The 
dashed curve shows calculated results using a constrain factor of 2.7 (plane 
strain) and the results are in excellent agreement with the experimental data, 
whereas, at the high-stress level, the predicted results gave shorter lives 
than the experimental data. At the high-stress level, the plastic-zone sizes 
are about a factor-of-4 larger than the low-stress level case, and plane-stress 
conditions may prevail. The dash-dot curve was calculated using an ct * 1.15. 
The calculated result'j are in better agreement with the experimental data than 
the results with an a •• 1.9. These results indicate that the constraint 
factor may vary with stress level and crack length. 

COMPARISON OF EXPERIMENTAL AND PREDICTED LIVES 
Figure 15 compares experimental (N^) and predicted (Np) lives for 
18 constant-amplitude load tests, 49 variable-amplitude load tests, and 
13 spectrum-load tests. The crack-closure model with a = 1.9 was used to 
predict crack-growth lives from the initial crack length c^^ to the final 
crack length c^. The ratio of predicted-to-experimental life (Np/Nj) ranged 
from 0.5 tc 1.8. The mean value of Np/Nj was 0.97 and the standard deviation 
was 0.27. 


CONCLUDING REMARKS 

An existing analytical crack-closure model (FAST) was used to correlate 
crack-growth rate data under constant-amplitude loading, and to predict crack 
growth under variable-amplitude and aircraft-spectrum loading. The model was 
based on the Dugdale model, but modified to leave plastically-deformed material 
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in the wake of the advancing crack tip. The leodel was used to calculate the 
crack-opening stresses as a function of crack length and load history under 
simulated plane-stress and plane-strain conditions. 

A previously developed crack-growth rate equation, in terms of Elber's 
effective stress-intensity factor range, threshold stress-intensity factor 
range, and fracture toughness, was used to correlate constant -amplitude rate 
data. The rate equation gives the "sigmoidal" shape connonly observed when 
fatigue crack-growth rate data is plotted against stress-intensity factor range. 
The five crack-growth constants in this equation were determined from constant- 
amplitude data on 2219-T851 aluminum alloy sheet material. The equation cor- 
related the constant-amplitude data over a wide range of stress ratios and 
stress levels quite well. 

The analytical closure model with a constraint factor of 1.9 and the rate 
equation were used to predict crack growth under variable-amplitude and 
aircraft-spectrum loading on the 2219-T851 aluminum alloy material. The proper 
constraint factor was determined from the constant-amplitude data. The model 
predicts the effects of load Interaction, such as retardation and acceleration. 
The ratio of predicted-to-experimental crack-growth lives (Np/N^) ranged 
from 0.5 to 1.8 in sixty-two variable-amplitude and spectrum load tests. The 
mean of Np/N^ was 0.97 and the standard deviation was 0.27. Thus, the 
analytical crack-closure model and the proposed crack-growth law predicted 
crack growth behavior in all tests quite well. 
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Spectrum for test M91 



FIG. 13 - Calculated crack-opening stresses as a function of crack length under 
typical aircraft spectrum loading. 


2.7 Plane strain 



u 


35 


FIG. U - Comparison of experimental [lA] and predicted crack-length-against-cycles 
curves for spectrum loading. 


test cycles 



FIG. 15 - Comparison of experimental [14] and predicted cycles to failure 
for 2219-T851 aluminum alloy material under constant-amplitude, 
variable-amplitude, and spectrum loading. 
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